Fluid Acceleration in the Bulk Turbulence of Dilute Polymer Solutions 
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We report experimental measurements of Lagrangian accelerations in the bulk of intense turbulent 
flows of dilute polymer solutions by following tracer particles with a high-speed optical tracking 
system. We observed a significant decrease in the acceleration variance in dilute polymer solutions. 
The shape of the normalized acceleration probability density functions, however, remain the same 
as in Newtonian water flows. We also observed an increase in the turbulent Lagrangian acceleration 
autocorrelation time in these dilute polymer solutions. 

PACS numbers: 47.27.Jv, 47.50.-d, 47.50.Ef, 47.57.Ng 



It has been known for many years that the addition 
of a very small amount of long-chain polymers to a fluid 
can have a dramatic effect on the turbulence near solid 
boundaries and result in significantly reduced flow resis- 
tance p. This "drag- reduction" phenomenon has been 
the subject of extensive studies 0, H, 0, 0, @, 0, S, 0- 
On the other hand, the effect of polymers on the bulk 
turbulence, away from boundaries, has received less at- 
tention. Tabor & De Gennes 0, Eo] and Lebedev & co- 



workers [Hi, [12, [13| proposed theories for homogeneous 
isotropic turbulence in dilute polymer solutions. Ex- 
periments 0, [H, 15 1 and numerical simulations 1^, 17 1 
have provided evidences that in bulk turbulence the small 
scales Eulerian quantities are damped by the polymer ad- 
ditives. Validation of these theories, however, has been 
very difficult due to the lack of detailed experimental 
data. 

Recently, there has been growing interest in studying 
turbulence from a Lagrangian viewpoint, i.e., by follow- 
ing the motion of the fluid particles [l8[. It has been 
shown, in particular, that the dynamics of turbulent flows 
with polymer additives is closely related to the correla- 
tion of velocity gradient following Lagrangian trajecto- 
ries [H. 

In this Letter, we report experimental measurements 
of fluid acceleration and its Lagrangian correlations in 
intense turbulent flows of dilute polymer solutions with 
various concentrations. The Taylor microscale Reynolds 
number of the turbulent flows before adding polymers 
was in the range of 140 < R\ < 485. We measured fluid 
accelerations by optically following tracer particles with 
a high-speed three-dimensional particle tracking system. 
We observed strong decreases of turbulent Lagrangian ac- 
celerations and an increase in acceleration anisotropy in 
these flows, similar effects on Eulerian small scale quanti- 
ties have been reported previously 14,1^. Moreover, we 



auto-correlation time of Lagrangian acceleration compo- 
nents increases with polymer concentration. These mea- 
surements provide detailed data of the effect of polymers 
on small scale Lagrangian turbulence. In addition, we ob- 
served an increase of the correlation time of Lagrangian 
acceleration magnitude, which suggests that larger scales 
of turbulence arc also affected. 

We studied a turbulent water flow generated between 
two counter-rotating coaxial baffled disks with blades 
that force the flow inertially rather than through a 
boundary layer, as shown in earlier experiments in a 
simflar apparatus M- The flow was seeded with 50/im 
or 26/im neutrally buoyant polystyrene particles. These 
particles have previously been shown to faithfully fol- 
low the turbulent flow in pure water at Reynolds num- 
bers up to 0(10'^) In this study, no difference was 
seen between measurements made with the two different 
sized particles in the water or in the polymer solutions. 
In addition, within the range of polymer concentrations 
studied in this work (0-lOppm), we observed only small 
changes in fluctuation velocity with polymer concentra- 
tion, comparable to the estimated uncertainties of veloc- 
ity measurement. 

A 35W pulsed frequency-doubled Nd:YAG laser was 
used to illuminate the tracer particles. A 4.1 x 4.1 x 2.05 
mm"^ volume was imaged onto four silicon strip detectors, 
each of which measured one coordinate of the particle 
position with spatial resolution 8/im/pixel. The silicon 
strip detectors can record up to 70,000 frames per sec- 
ond. The recording frequencies were chosen such that the 
time resolution of the measurement is at least 100 frames 
per Tri, sufRcient for accurate acceleration measurement. 
The details of the flow and imaging system have been 
described in 
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observed that both effects increase with polymer concen- 
tration. The probability density function (PDF) of accel- 
eration normalized by the standard deviation, however, 
is not affected by the polymer additives. The measured 



The measured particle positions were connected to 
form three dimensional particle trajectories using the al- 
gorithms described in Ref. [2lj. Velocities and acceler- 
ations were computed by convolution of the trajectories 
with an appropriate kernel based on successive deriva- 
tives of a Guassian filter of width ak, which performs 
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simultaneously a low-pass filtering and the differentia- 
tion [2^. The variances of acceleration components, (a,^), 
for both water and the polymer solutions, were estimated 
(of) is computed as a function of ak and 



as m 
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I.e., 



then is extrapolated to ak going to zero, in order to 
remove the effect of measurement noise on acceleration 
variance. 

We developed a protocol to prepare the polymer so- 
lutions such that the experiments are robust and repro- 
ducible. First a high concentration stock solution was 
prepared by mixing the dry polymer powder with wa- 
ter. A magnetic stir was used to mix the solution at 
slow speed for 16 hours. Then the stock solution stood 
for several hours and was depressurized to remove solu- 
ble gases. Finally it was gently added to the turbulence 
chamber to form a dilute solution at the desired con- 
centration. We typically waited 10 to 15 minutes before 
starting the measurement so that the solution was well- 
mixed. In this Letter, we present measurements on solu- 
tions of polyacrylamide (PAM) 18522 from Polysciences. 
All initial stock solutions were 400 ppm except one ex- 
periment with a final concentration of Ippm, which had 
a stock solution of only 20 ppm. 

According to the manufacturer, the weight average 
molecular weight, M^,, for the 18522 PAM is M„ = 
18 X lO^amu. The manufacturer also gave a measure of 
the polydispersity of the sample M^/Mn = 3 — 4, where 
Mn is the number average molecular weight. Using the 
Zimm model Q we calculated that the radius of gyra- 



tion Rg = 0.5/Ltm, the maximum extension 
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and the polymer relaxation time Tp = 43ms for this poly- 
mer. As it is widely accepted, the most important dimcn- 
sionless number that characterizes the polymer dynamics 
in turbulence is the Weissenberg number Wi = Tp/Tr,, 
which measures the ratio of the polymer relaxation time 
to the Kolmogorov time r^, the smallest time scale of 
turbulence. Only when Wi > 1, the polymer molecules 
can be stretched by the turbulent flow and may affect the 
turbulence. The corresponding Weissenberg number for 
each Reynolds number that we have measured is shown 
in Table [D Some experiments have shown that the effect 
of polymer on turbulence is sensitive to the polydisper- 
sity of the polymer sample with the longer polymers in 
the sample having a greater effect on the turbulence [2^ . 
In view of this, the Tp and Wi calculated from the M^, 
may be treated as a lower bound. 

Figure [1] shows the component acceleration PDF for 
water and the polymer solution of 5ppm by weight at 
R\ ~ 285. The unnormalizcd PDFs for the dilute poly- 
mer solution are strongly suppressed, which indicates 
damping of dissipation scale fluctuations by the presence 
of polymers. It is interesting to note that the PDF 's of 
the normalized acceleration, ai/(a|)^/^, for water and for 
the dilute polymer solution overlap with each other, as 
shown in Figures [T] (c) and (d), which suggests that the 
accelerations arc suppressed only by a scaling factor. 
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TABLE I: Parameters of the experiments, fp is the rota- 
tion speed of the baffled disks. Rx is the Taylor microscale 
Reynolds number. 77 and are the Kolmogorov length and 
time scales, respectively. R\, 77, and r,, refer to the turbulent 
water flows before adding polymer solutions, and are mea- 
sured as in Ref. [2l|]. Nf is the temporal resolution of the 
recording system, in frames per r,,. Wi = Tp/rn is the Weis- 
senberg number. 



The rms acceleration, (af)^^^, in the polymer solutions 
is dramatically lower than in water. Figure [2] shows the 
ratio of the rms acceleration in polymer solutions to that 
of water for the same rotation rate of the propellers. In 
the range of polymer concentrations investigated, the de- 
pendence on concentration is roughly linear. There is no 
clear dependence on R\ or Wi. 

It should be pointed out that the decrease in accelera- 
tion variance can not be explained by a decrease in energy 
dissipation rate e. In earlier experiments carried out in 
similar geometry, it has been shown that the total energy 
dissipation rate remains unchanged when polymers were 
injected into the flow if baffled disks were used to pro- 
vide the driving (20| . The incrtial forcing by the vanes 
was not altered by the presence of polymers. 

The anisotropy of the small scales is also increasing 
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FIG. 1: (color online). PDF of acceleration components at 
R\ = 285. o are for water. • are for a 5ppm polymer so- 
lution. The width of the filter used to calculate acceleration 
is cTfc = 0.23r^. (a) P(a,), (b) P{a^), (c) P{a^/(aiy^^) (d) 
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FIG. 2: (color online). Ratio of the rms acceleration in poly- 
mer solution to that of water vs. concentration of the poly- 
mer solution. □, Rx — 140, axial component; o, = 200, 
axial component; •, R\ = 200, horizontal component; A, 
Rx = 285, axial component; A, Rx ~ 285, horizontal compo- 
nent; o, Rx ~ 485, axial component; For the data point at 
Rx = 485 and concentration 5ppm, statistics from the first 3.3 
hours of the run were used. The variance changes less than 
10% over this time. Data points at fp — 0.6 hz and concentra- 
tion 5 ppm were shifted to 5.25 ppm for clarity. INSET: Ratio 
of the horizontal acceleration variance to the axial accelera- 
tion variance as a function of concentration of the polymer 
solution. •, Rx = 200; A, i?A = 285. 



upon addition of the polymer. The insert in Figure 
shows the ratio of the horizontal acceleration variance to 
the axial acceleration variance as a function of polymer 
concentration. If the polymers were causing the energy 
cascade to be truncated at a smaller wave number, then 
the separation between large and small scales would be 
reduced and the anisotropy of the large scales would af- 
fect the small scales more. This is similar to seeing the 
anisotropy increase with decreasing R\. The increase 
of anisotropy might also come from the alignment of 
polymers in particular directions due to the large scale 
anisotropic flow. Since the polymer responses in longitu- 
dinal and transverse directions (relative to the alignment 
of polymers) are different [2^, this difference could ap- 
pear as anisotropy in acceleration. 

Since the dynamics of polymer solutions depends 
on the velocity field following Lagrangian trajectories, 
we studied the Lagrangian acceleration autocovariance 
Cy (r) = {ai(t)aj{t + t)) and the Lagrangian autocorre- 



lation Cij{T) 



{a,{t)aj{t+r}) 



To account for the 



effect of the finite volume bias, we calculated the modi- 
fied autocovariance Cij (r) 



{a^{t)a^{t + t)) 



t&CT 



Htf),, (1) 



in which Of&r indicates that only measurements of accel- 
eration where a measurement at both t and t + r existed 



on a particle track were averaged over and ()t indicates 
that all measurements of acceleration were averaged over. 
The autocorrelation is then 



G^J{T) 



(a|)V2 



(2) 



where {of) is determined from the procedure described 
above. 

Figure [3] shows Cij (r) and (r) for water and poly- 
mer. The height of the autocovariance is depressed 
because, as we have discussed above, (a|) is dramati- 
cally decreased in the polymer solutions. It has been 
shown that the integral of the acceleration correlation 
vanishes [i^l and the zero-crossing time of the correla- 
tion curve is usually used as a measure of the correlation 
time [1^ [2^ . The curve at the zero-crossing has a shal- 
low slope. Small measurement errors in the acceleration 
correlation can cause large deviations in the zero-crossing 
time. A more robust measure is the exponential decay 
time Tx/e, defined such that C(T]^/e) — V*^- pure wa- 
ter, the Ti/e times are approximately 0.7-0.8 times the 
Kolmogorov times. As shown in Figure [31 in the range 
of polymer concentration used, the ti/^ for polymer is 
about 3-4 times the ti/^ for water and increases with 
polymer concentration. 

The width of the acceleration correlation is propor- 
tional to the Kolmogorov time, r^. In pure water we ob- 
serve that component acceleration correlation functions 
for R\ = 285 to 970 collapse when the time is normal- 
ized by T,j (30! . The Kolmogorov time in turn is propor- 
tional to v^/'^ . The acceleration variance is proportional 
to i'"^/^. If the polymer solutions were to behave simply 
as if they had an increased effective viscosity, then we 
would expect the decrease in the acceleration variance 
to match the increase in the correlation time. Instead 
we find that the decrease in (a|) is much larger than 
the increase in the autocorrelation time. The accelera- 
tion variance for a 3.5 ppm solution is more than 4 times 
higher than for a 10 ppm solution. However the Ti/g for 
the 10 ppm solutions is not even twice as much as that 
of the 3.5 ppm solution. The backreaction of the poly- 
mer is obviously more complex than just a change in the 
bulk viscosity and is most likely anisotropic and depen- 
dent on the local velocity gradients [l^ as for example 
the elongational viscosity is known to change by orders 
of magnitude but not the shear viscosity 0. 

We also measured the autocorrelation of the mag- 
nitude of Lagrangian acceleration, defined as C|a| — 
( ( I ~ ^ I "^1 j ) ( I ^j]^^.}}r' ^ I "^1 '' ^ ^ and obtained using an estima- 
tor equivalent to Eq. ([2]) . The results are shown in Fig. |4l 
The time scale given by the acceleration magnitude cor- 
relation increases with polymer concentration. This may 
indicate that the large scales of the flow are affected by 
the polymer as well, since the acceleration magnitude is 
correlated over much longer time scales of the flow than 
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FIG. 3: (color online). Autocovariance of the axial acceler- 
ation component for R\ = 285 at different polymer concen- 
trations. • water; ▲ lOppm; A 5ppm; □ 3.5ppm. INSET: 
Corresponding autocorrelations. 



the acceleration components. Due to the finite size of the 
measurement volume, we only measured the correlation 
of the acceleration magnitude over a few r,,. ft is possible 
that, at longer times, the curve for the polymer solution 
rejoins that for water. 

In summary, we studied the effect of polymers on the 
small scales of turbulence by measuring the fluid accel- 
eration in turbulent flows of dilute polymer solutions. 
We observed strong decreases of acceleration variance 
in polymer solutions while the shape of the acceleration 
PDF remains the same as in water. The effect on accel- 
eration variance was insensitive to the Reynolds number 
or the Weissenberg number. The Lagrangian correlation 




FIG. 4: (color online). C\a\ (j) for water and polymer, o water 
CTfc = 0.35r,,. • water R\ = 285, = 0.31r,,. □ 



i?A = 485, 
5ppm, Rx 



485, (Tfe = 0.35t^. a 5ppm, R\ 



285, 



Ok 



0.31r^. ▲ lOppm, 7?a = 285, gu = 0.31r„. The polymer 
statistics for R\ = 485 were taken from the first 3.3 hours of 
the run. 



time of acceleration components increases with polymer 
concentration. The observed effect of polymer additives 
cannot be explained by a change of viscosity alone. In 
addition, the measured Lagrangian correlation time of ac- 
celeration magnitude increases with concentration, which 
implies that polymers may affect larger scales. 
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